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Abstract 

The spin-hah Heisenberg antiferromagnet (HAF) on the square and triangu- 
lar lattices is studied using the coupled cluster method (CCM) technique of 
quantum many-body theory. The phase relations between different expansion 
coefhcients of the ground-state wave function in an Ising basis for the square 
lattice HAF is exactly known via the Marshall-Peierls sign rule, although no 
equivalent sign rule has yet been obtained for the triangular lattice HAF. 
Here the CCM is used to give accurate estimates for the Ising-expansion co- 
efhcients for these systems, and CCM results are noted to be fully consistent 
with the Marshall-Peierls sign rule for the square lattice case. For the trian- 
gular lattice HAF, a heuristic rule is presented which fits our CCM results for 
the Ising-expansion coefficients of states which correspond to two-body exci- 
tations with respect to the reference state. It is also seen that Ising-expansion 
coefficients which describe localised, m-body excitations with respect to the 
reference state are found to be highly converged, and from this result we infer 
that the nodal surface of the triangular lattice HAF is being accurately mod- 



eled. Using these results, we are able to make suggestions regarding possible 
extensions of existing quantum Monte Carlo simulations for the triangular 
lattice HAF. 
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I. INTRODUCTION 



The coupled-cluster method [Q-^ (CCM) has been previously applied 
||]"0|-[2^ to a number of unfrustrated and frustrated lattice quantum spin sys- 
tems with a great deal of success. Recently, it has been shown for a CCM 
calculation, within the so-called SUB2 approximation which contains all two- 
body correlations, that the ket-state correlation coefficients can be related to 
the Marshall-Peierls sign rule for the square lattice Heisenberg antifer- 



romagnet (HAF). Furthermore, it was shown in this treatment |24| that the 
CCM ket-state correlation coefficients for the Heisenberg model on the trian- 
gular lattice also demonstrate a seemingly regular pattern. This relationship 
between the Marshall-Peierls sign rule and the CCM ket-state correlation coef- 



ficients was further investigated [25| and clarified for a spin model with both 
nearest-neighbour and next-nearest-neighbour exchange, namely the J1-J2 
model, on the linear chain and square lattices. The Marshall-Peierls sign rule 
at the Heisenberg point on the square lattice was found to be preserved by all 
of the Ising-expansion coefficients within a localised approximation scheme, 
as well as within the SUB2 scheme. No assumption was made other than the 
choice of an initial reference state with respect to which the ground-state wave 
function was approximately constructed in the infinite lattice limit using the 
CCM formalism. It was also possible to obtain a quantitative value for the 
point at which the sign-rule breaks down with increasing antiferromagnetic 
next-nearest-neighbour exchange for the square lattice case. This quantitative 
value was found to be in good agreement with exact diagonalisation calcula- 



tions |27|. 

The Marshall-Peierls sign rule for Heisenberg antiferromagnets on bipar- 
tite lattices is an exact statement for the phase relations between the expan- 
sion coefficients of the ground-state wave function in an Ising basis. This 
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statement is of intrinsic interest because it provides exact information regard- 
ing the ground-state wave function in cases (e.g., the square lattice HAF) for 
which no general exact solution has yet been determined. It is also of interest 
because by using this rule one implicitly has full knowledge of the 'nodal sur- 
face' of the ground-state wave function within this basis. (Note that whenever 
the words 'nodal surface' are referred to in this article, this is taken to mean 
the magnitudes and signs of these expansion coefficients.) Of course, for the 
square lattice HAF the signs are exactly known via the Marshall-Peierls sign 
rule, although the magnitudes are known only via exact or approximate cal- 
culation. Quantum Monte Carlo (QMC) calculations for fermionic and spin 
systems demonstrate the infamous "sign-problem", and so such information 
regarding the signs of the expansion coefficients is of importance in these 
calculations. 

For non-bipartite lattices (e.g., the triangular lattice HAF) or for other 
models containing frustration (e.g., the J1-J2 model on the square lattice) 
corresponding exact sign-rule theorems are generally not known. Hence, var- 
ious approximate ways of simulating the nodal surface have been developed 
for use, for example, in fixed-node quantum Monte Carlo calculations. How- 
ever, one is always limited in such simulations by not knowing how close is 
the trial wave function or guiding wave function that one uses to the true 
ground-state wave function. It is therefore of considerable interest to develop 
theories which simulate the nodal surface from an ab initio standpoint, such as 
is done here via the CCM. In this way, one might be able to utilise this CCM 
information concerning the nodal surface in a QMC calculation, or at least 
suggest possible extensions of previous QMC calculations. It is also conceiv- 
able that one might observe patterns in the expansion coefficients obtained 
via the CCM and so infer rules concerning these coefficients or sub-classes of 
them which relate to specific types of excitation. 
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In this article, we briefly introduce the CCM formahsm for the lattice 
quantum spin systems and also make an explicit link between the expansion 
coefficients of the ground-state wave function and the CCM ket-state corre- 
lation coefficients. This link is investigated for the square and triangle lattice 
HAFs, and the Ising-expansion coefRcients corresponding to states which cor- 
respond to two-body excitations with respect to a reference state are deter- 
mined within well-defined approximation schemes. These coefficients are seen 
to be well-converged using a localised approximation scheme. For the trian- 
gular lattice HAF, a heuristic rule is presented which fits our CCM results for 
the Ising-expansion coefficients of states which correspond to two-body exci- 
tations with respect to the reference state. The Ising-expansion coefRcients 
of states which correspond to m-body excitations are also shown, within this 
same localised approximation scheme, to be well-converged. Patterns are also 
seen in the Ising-expansion coefficients which correspond to these m-body ex- 
citations for the triangular lattice HAF, although no generic sign rule for them 
is formulated. Finally, suggestions are made regarding possible extensions of 
existing quantum Monte Carlo simulations for the triangular lattice HAF. 

II. THE SPIN MODEL AND THE CCM FORMALISM 

A. The Heisenberg Model 

In this paper we consider the zero-temperature (T = 0) properties of the 
spin-half Heisenberg antiferromagnet (HAF) quantum spin system, which is 
described by the Hamiltonian 



where the index i runs over all N lattice sites on the square and triangular 
lattices with periodic boundary conditions, and the index j runs over all 
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(1) 



nearest-neighbour sites to i. The angular brackets indicate that each nearest- 
neighbour bond (or link) is counted once and once only. 

The Heisenberg model on the square lattice has not been solved exactly 
up till now, although it been extensively studied using various approximate 
methods [p^31|. Runge Q| has performed the most accurate Monte Carlo 



simulation to date for the square-lattice, isotropic HAF. He finds a value for 
the ground-state energy per spin of £'g/A^=— 0.66934(4), and a value for the 
sublattice magnetisation which is 61.5%ib0.5% of the classical value. Exten- 



sive CCM calculations have also been carried out for this model [23[ giving a 
value for the ground-state energy per spin of £'g/A^=— 0.66968, and a value for 
the sublattice magnetisation which is 62% of the classical value. In compari- 
son, linear spin-wave theory (LSWT) gives a value of Eg/N=— 0.658 for 
the ground-state energy, and a value for the sublattice magnetisation which 
is 60.6% of the classical value. 

Again, no exact solution exists for the Heisenberg model on the tri- 
angular lattice, although many approximate calculations have been carried 
out [ 23, 3^ 37 1 . Exact series expansion |33| calculations obtain a value for 
the ground-state energy per spin of this model of Eg/N=— 0.551. Simi- 
larly, spin-wave theory (SWT) ||3^ gives a value of Eg/N=— 0.5388, and a 
fixed-node quantum Monte Carlo (FNQMC) |^5| calculation gives a value of 
Eg/N=— 0.5431 lb 0.0001. Exact diagonalisations of finite-sized clusters of 
spins |^6| which have been extrapolated to the infinite lattice limit give a 
value for the ground-state energy per spin of Eg/N=— 0.5445. In comparison, 
a recent coupled cluster method (CCM) calculation |2^ predicts a value for 
the energy of Eg/N=— 0.5505, which is fully consistent which the best of the 
other results obtained for this system so far. 

Classically, this system orders as a Neel state on the triangular lattice 
such that each nearest-neighbour pair of spins makes an angle of 120° to each 

6 



other. The results of approximate theories for the amount of this ordering 
that remains in the quantum hmit, i.e., the sublattice magnetisation, are 



typicahy |23,p4,p6| that 50% of the classical ordering remains in the quantum 



limit. The FNQMC calculation of Ref. [35| places this value at 60%, which 
is probably still too high. A notable exception to all of these results quoted 



above is the result of exact series expansion calculations [33| which predicts 
that as little as 20% of the classical ordering remains. 

B. The CCM Formalism 



The CCM formalism [yj-^ is now briefly considered, although the reader 
should note that more detailed descriptions of the CCM applied to spin sys- 
tems are given in Refs. [pX|-25|. To calculate the ground state wave function 



l^*) of a spin system we start with a model state \ ^) and a correlation operator 
S such that 

|^) = e^|$) , (2) 

where the ket-state correlation operator S may be written as, 

S = Y.SjCt . (3) 

The correlation operator S is formed from a linear combination of multi- 
configurational creation operators, {C^}, (which themselves are formed from 
products of spin raising operators) multiplied with the relevant ket-state cor- 
relation coefficients {Si}. In addition, we define a set of destruction operators, 
{Cy}, which are the Hermitian adjoints of {Cf}. 

For the square lattice case, the model state is chosen to be the classical 
Neel state. The lattice is divided into two sublattices and one sublattice is 
populated with 'up' spins and the other with 'down' spins. In order to treat 
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the spins equivalently the local spin axes of the 'up' spins are rotated by 180° 
about the y-axis, which is mathematically written as, 



-s^ ; sy^sy ; 



The Hamiltonian may now be written, with the introduction of an anisotropy 
coefficient x on the off-diagonal elements of H, as 

(id) 

where x = 1 corresponds to the isotropic Heisenberg Hamiltonian of Eq. 
dH). The Marshall-Peierls sign rule for the original Hamiltonian of Eq. (|l|), 
in terms of rotated spin coordinates, now differs from the original sign rule 
|26|. For the ground-state wave function, j'l') = J2i ^/ 1-^)) the Ising-expansion 



coefficients {^/} must now all be greater than (or equal to) zero with respect 
to a complete set of Ising basis states {|/)} in the rotated spin coordinates. 

For the triangular lattice case, we choose \^) to again be the classical 
Neel state, and we start the CCM calculation by dividing the lattice into 
three sublattices, denoted {A,B,C}. The spins on sublattice A are oriented 
along the negative z-axis, and spins on sublattices B and C are oriented at 
+120° and —120°, respectively, with respect to the spins on sublattice A. 
In order both to facilitate the extension of the isotropic HAF to include an 



Ising-like anisotropy first introduced by Singh and Huse [33| and to make 
a suitable choice of the CCM model state, we perform the following spin- 
rotation transformations. Specifically, we leave the spin axes on sublattice A 
unchanged, and we rotate about the y-axis the spin axes on sublattices B and 
C by -120° and +120° respectively, 

^ X ^ ^ I ^ V3 , 

^^B ^ ^^B ' ^c^^c ' 

Z X ^ z z X ^ z 



We may rewrite Eq. (||) in terms of spins defined in these local quantisation 
axes for the triangular lattice with a further introduction of an anisotropy 
parameter A for the non-Ising-like pieces, 

where A = 1 corresponds to the isotropic Heisenberg Hamiltonian of Eq. ([l|). 
We note that the summation in Eq. (^) again runs over nearest-neighbour 
bonds, but now also with a directionality indicated by {i j), which goes 
from A to B, B to C, and C to A. We note that no exact sign rule has yet 
been proven for the Hamiltonian of Eq. (^). When A = 0, the Hamiltonian in 
Eq. (0) describes the usual classical Ising system with a unique ground state 
which is simply the fully aligned ("ferromagnetic") configuration in the local 
spin coordinates described above. 

From the Schrodinger equation, H\^) = E\^), we obtain an expression 
for the ground state energy which is given by 

This equation illustrates one of the key aspects of CCM, which is namely the 
similarity transform - here of the Hamiltonian. The similarity transform of 
any operator may be expanded as a series of nested commutators such that 
an exact expression for the ground-state energy in terms of the ket-state cor- 
relation coefficients, {Sj}, which is already determined in the infinite lattice 
limit ^ oo, can often be found for whatever (non-trivial) approximations 
are made for S. The approximation in S means that the we obtain approx- 
imate values for ket-state coefficients, {5/}, and so the ground-state energy 
is approximately determined, although, as stated above, the equation for the 
ground-state energy is an exact expression of one or more of these coefficients. 
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To determine the CCM ket-state coefficients we operate on the Schrodinger 
equation with exp(— S") and then by {^{Cy , for a given cluster configuration 
indicated by the index /. An explanation of how the CCM equations may be 
derived and then solved is presented in Refs. 

The three most commonly employed approximation schemes are: (1) the 
SUBn scheme, in which all correlations involving only n or fewer spins are 
retained, but no further restriction is made concerning their spatial separa- 
tion on the lattice; (2) the SUBn-m sub-approximation, in which all SUBn 
correlations spanning a range of no more than m adjacent lattice sites are 
retained; and (3) the localised LSUBm scheme, which retains all multi-spin 
correlations over distinct locales on the lattice defined by m or fewer contigu- 
ous sites. Note that all 30 cluster configurations which are used in the LSUB4 
approximation for the triangular lattice antiferromagnet, and which are inde- 
pendent under the six-point symmetry group inherent in the Hamiltonian of 
Eq. (0), are hsted in Fig. |l|. 

We now consider an expansion of the ground-state wave function in a 
complete Ising basis {|/)} (in terms of the local coordinates after rotation). 
This may be again written as, |^') = J2i^i\^)i where the sums over / goes 
over all 2^ Ising states, and we find that this expression naturally leads from 
Eq. (^) (also see Ref. [^) to an exact mapping of the CCM correlation 
coefficients {Sj} to the Ising-expansion coefficients {^'/}, which is given by 

It is possible to match the terms in the exponential to the 'target' configura- 
tion of CY in Eq. (|9|), and so obtain a numerical value for the {^'/} coefficients 
once the CCM ket-state equations have been derived and solved for a given 
value of the anisotropy. As the target configuration is formed from a finite 
number of spin lowering operators, there is a cut-off point in the Taylors se- 
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ries expansion of the exponential above which no contribution to ^/ in Eq. 
(|9|) can possibly occur. This matching may be achieved in two ways: 
analytically at the SUB2 level of approximation or at small LSUBm levels 
of approximation; or by using computational-algebraic techniques. We note 
that we use, for the sake of consistency, the same configurations for the 'tar- 
get' configurations in in Eq. (^) as are used in the ket-state correlation 
operator S for the approximation schemes define above. We also note that, 
although this is an approximate calculation, we are already dealing with the 
infinite lattice limit, N ^ oo. 

The next section describes our results for the set of {^i} coefficients for 
the spin-half Heisenberg model on the square and triangular lattices using the 
CCM formalism, and discusses possible rules that can be inferred from these 
results. 

III. RESULTS 

A. The Two-body Excitations 

The results for the Ising-expansion coefficients {^i} corresponding to two- 
body excitations with respect to the model state for the spin-half square 
lattice HAF {x = 1) are shown in Fig. ||. (Note that the solution to the CCM 
equations is tracked from x = 0, at which point we know that all of the ket- 
state correlation coefficients {Si} are zero, to the isotropic Heisenberg point.) 
The results for the two-body {^i} coefficients are then determined and it is 
found that these coefficients are all positive for the LSUBm approximation 
scheme (with m < 8) and within the SUB2 approximation. (This result for 
the SUB2 case was previously seen in Ref. [^.) This is an expression of the 
Marshall-Peierls sign rule ||2^ which provides an exact relation between the 
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{^i} expansion coefficients. However, we note that the CCM predicts such 
behaviour with no recourse to exact proofs, and this prediction is a natural 
result of the CCM calculations. Furthermore, the convergence of the LSUBm 
series of results compared to SUB2 calculation can clearly be seen in Fig. ^. 
For short lattice distance, we believe that the LSUBm results are much better 
model of the 'nodal surface' of states corresponding to two-body excitations 
than those of the SUB2 calculation. 

For the spin-half triangular lattice HAF (A = 1), the results for the Ising- 
expansion coefficients {^'/} corresponding to two-body excitations with re- 
spect to the model are shown in Fig. ^ and in Table 1. (Note again that the 
solution to the CCM equations is tracked from A = 0, at which point again 
we know that all of the ket-state correlation coefficients {Sj} are zero, to the 
isotropic Heisenberg point.) A rule that fits the behaviour seen in Fig. ^ can 
now be stated as: the coefficients are found to be positive if the lattice vec- 
tor of between the two spin-raising operators in S connects sites on different 
sublattices; conversely, the {^/} coefficients are found to be negative if the 
lattice vector connects sites on the same sublattice. It is again seen that the 
LSUBm series is clearly well-converged for small lattice distance R. This se- 
ries of LSUBm results also clearly present a much better representation of the 
'nodal surface' for the two-body excitations, for small lattice distance than 
the corresponding {^/} coefficients within the SUB2 approximation [p^]. 



B. The m-body Excitations 

For the Ising-expansion coefficients {^'/} corresponding to m-body exci- 
tations for the square lattice HAF (x = 1) , it has previously been noticed |25| 



that the CCM calculation outlined above predicts that they are all positive. 
This is an indication that the CCM results are (once more) fully consistent 
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with the exact Marshall-Peierls sign rule for this model. It should be again 
noted that the only assumption made is the initial model state, and that the 
rest follows naturally on from the CCM calculation. 

Table 1 indicates the {^i} expansion coefficients for for the triangu- 
lar lattice HAF (A = 1) states corresponding to localised, m-body (with 
m < 4) configurations (i.e., those also used within the LSUB4 approxima- 
tion) with respect to the model state. One should note however that these 
{^i} expansion coefficients are determined within the SUB2 and LSUBm 
(with m = {3,4,5,6}) approximations. It is again seen in Table 1 that all 
of the expansion coefficients are well converged for the LSUB5 and LSUB6 
levels of approximation. Note, however, that within the LSUB5 and LSUB6 
approximation many other {^i} coefficients have been determined, within 
the 'consistency' assumption explained above. 

Patterns in the {^i} expansion coefficients for the triangular lattice HAF 
(A = 1) for states corresponding to m-body excitations (with m > 1) have 
also been observed. It was seen that, with m odd, all of the coefficients 
for excitation clusters which are equivalent under both the symmetries of 
the lattice (a point group of order twelve) and the Hamiltonian (a point 
group of order six), were found to have both CCM correlation coefficients 
and Ising-expansion coefficients which are zero at all levels of approximation; 
e.g., configurations (1), (3), and (5) in Fig. |l] and in Table 1. In contrast, 
those {^/} coefficients corresponding to configurations with odd numbers of 
spins which are equivalent under the symmetries of the lattice but are not 
equivalent under the symmetries of the Hamiltonian were found to have equal 
magnitudes but opposite signs. That is, there exists pairs (12/6 = 2) of 
coefficients, for the odd-number spin-excitation clusters, which differ with 
respect to each other only by a phase factor of —1; e.g., those pairs (7,12), 
(8,10), (16,19), (18,25), (21,23), and (27,29) in Fig. g and in Table 1. In 
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contrast, for the Ising-expansion coefficients {^i} corresponding to m-body 
excitations, with m even, we find that those {^i} coefficients corresponding 
to configurations which are equivalent under symmetries of the lattice but are 
not equivalent under the symmetries of the Hamiltonian are exactly the same. 
That is, there exists pairs (12/6 = 2) of coefficients, for the even-number spin- 
excitation clusters, which have a phase factor with respect to each other of 
unity; e.g., those pairs (9,13), (22,26), and (28,30) in Fig. | and in Table 1. 

We note that these simple relations were obtained with no assumption 
other than the model state, which is of course the classical ground state of 
this system. This is in contrast to variational calculations spin wave 



theory |3^ and FNQMC |35] which assume that the classical phase relations 



are correct phase relations for their quantum mechanical counterparts. It 



has been shown |18] that in a case in which the Hamiltonian shares some 
symmetry of the model state then the values of the CCM correlation coeffi- 
cients in the ground state, for configurations which are equivalent under these 
symmetries, must be completely identical. We believe that the 'rules' pre- 
sented above for the m-body expansion coefficients for the triangular lattice 
HAF are therefore a refiection of the symmetries inherent in the Hamiltonian 
and lattice. Furthermore, it remains an open question as to the solution of 
all branches of coupled, non-linear, ket-state CCM equations for high-order 
calculations. One should note that it might be possible to obtain another so- 
lution branch which might not demonstrate the full symmetries of the lattice 
and Hamiltonian, as has been seen for other systems treated by CCM before. 
A full treatment of this subject would form the contents of another article. 

Apart from these simple observations, we cannot infer a more general rule 
as to the behaviour of the expansion coefficients for m-body excitations with 
m > 2. However, we do have results for the signs of all of the expansion co- 
efficients within the LSUBm approximation scheme, as is illustrated in Table 

14 



1, and hence the CCM is simulating, from an a priori view-point, the nodal 
surface of this model. From the amount of convergence of the results of the 
expansion coefficients within the LSUBm scheme described above, we there- 
fore also believe that we have an accurate simulation of the nodal surface of 
this model using the LSUBm approximation scheme. 



C. Relevance to QMC Calculations 

For the triangular lattice HAF, a 'naive' picture of how important each 
excitation configuration can be found by looking at the magnitude of its cor- 
responding coefficient, illustrated in Table 1. In this manner, we find that 
the two-body, nearest-neighbour configuration (configuration (2) in Fig. |l|) 
is the most important correlation for this model. By this rationale, the next 
most important configurations are a 'diamond-shaped', four-body configura- 
tion (configuration (11) in Fig. [l|) and a 'dog-leg', three-body configuration 
(configurations (8) and (10) in Fig. [l|). A previous FNQMC calculation |35| 



contains all two-body correlations and this particular three-body correlation, 
and so a straightforward extension of this calculation would be to include 
this four-body configuration. However, from this 'naive' interpretation of our 
results it appears that, as well as the two-body configurations and the two con- 
figurations mentioned above, many other localised three-body and four-body 
configurations are also very important. (These configurations are localised 
in the same sense of those configurations contained in the LSUBm approx- 
imation.) Hence, to attempt an accurate simulation of the this model one 
probably needs to find a way of also including many, if not all, of these lo- 
calised, higher-order correlations. One way of achieving this would therefore 
be to use the nodal surface predicted by the CCM in a FNQMC calculation. 
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IV. CONCLUSIONS 



In this article it has been shown that the expansion coefficients of the 
ground-state wave function in the infinite lattice limit may be approximately 
determined via CCM calculations. The two-body expansion coefficients have 
thereby been obtained as a function of lattice distance. For the triangular 
lattice HAF, a heuristic rule which fits CCM results for the coefficients cor- 
responding to two-body excitations with respect to the model state was also 
obtained. Certain patterns, for the triangular lattice HAF, were also seen for 
coefficients corresponding to higher-order excitations, although no complete 
heuristic rule for all Ising states could be inferred. For the square lattice HAF, 
it has once more been noted that all of these expansion coefficients obey the 
Marshall-Peierls sign rule for the LSUBm and SUB2 approximations with no 
other assumption other than the model state. The convergence of the Ising- 
expansion coefficients for the LSUBm series of results was also seen for small 
lattice distance for both lattices. In this manner, the LSUBm results are 
believed to provide an accurate simulation of the true nodal surface of these 
models at small lattice separation. 

Finally, suggestions were made with respect to further QMC calculations 
for the triangular lattice HAF. The first was the inclusion of a 'diamond- 
shaped' cluster in a FNQMC calculation, which we believe to be the most 
important cluster not yet used in such a calculation. However, we also infer 
from our results that many three- and four-body clusters are important in 
order to accurately simulate this model. One method of including these cor- 
relations in a QMC calculation would be to use the nodal surface predicted 
by CCM calculations. 
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FIG. 1. Excitation cluster configurations for the triangular lattice antifcrromagnet within the 
LSUB4 approximation. Each hexagon marks the lattice position of a spin raising operator applied 
to the model state at a given point on the triangular lattice. 
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FIG. 2. Results for the Ising-expansion coefficients, plotted as a function of the lattice distance 
R, corresponding to two-body excitations with respect to the model state for the spin-half, square 
lattice HAF {x = 1) obtained via the LSUBm approximation scheme (with m = {4, 6, 8}) and the 
SUB2 approximation. 
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FIG. 3. Results for the Ising-expansion coefficients, plotted as a function of the lattice distance 
R, corresponding to two-body excitations with respect to the model state for the spin-half, trian- 
gular lattice HAF (A = 1) obtained via the LSUBm approximation scheme (with m = {3, 4, 5, 6}) 
and the SUB2 approximation. 



23 



Table 1. Expansion coefficients {^i} for the spin-half, triangular lattice HAF 
(A = 1) for the clusters given in Fig. |l] within LSUBm (with m = {3, 4, 5, 6}) 
approximation scheme and the SUB2 approximation. For the sake of consis- 
tency, the only clusters expansion coefficients {^i} that have been determined 
are those which are also used within the given approximation. 



Cluster 


SUB2 


LSUB3 


LSUB4 


LSUB5 


LSUB6 


1 


0.000000 


0.0000000 


0.000000 


0.000000 


0.000000 


2 


0.114444 


0.1280959 


0.141570 


0.145502 


0.149245 


3 


- 


0.0000000 


0.000000 


0.000000 


0.000000 


4 


0.000941 


0.0026416 


0.007384 


0.007310 


0.008302 


5 


- 


0.0000000 


0.000000 


0.000000 


0.000000 


6 


-0.008838 


-0.024744 


-0.028535 


-0.031743 


-0.032843 


7 


- 


- 


0.004098 


0.004380 


0.004694 


8 


- 


-0.029600 


-0.039855 


-0.044560 


-0.047332 


9 


- 


- 


0.020290 


0.021182 


0.022933 


10 


- 


0.029600 


0.039855 


0.044560 


0.047332 


11 


- 


- 


0.055163 


0.060557 


0.065889 


12 


- 


- 


-0.004098 


-0.004380 


-0.004694 


13 


- 


- 


0.020290 


0.021182 


0.022933 


14 


- 


- 


0.014343 


0.013776 


0.014106 


15 


0.000428 


- 


-0.001568 


-0.002414 


-0.002185 


16 


- 


- 


0.000133 


0.000437 


0.000409 


17 


0.004417 


- 


0.007473 


0.008187 


0.009555 


18 


- 


- 


-0.001828 


-0.002663 


-0.002865 


19 






-0.000133 


-0.000437 


-0.000409 


20 






0.019778 


0.020721 


0.022046 


21 






0.005219 


0.006398 


0.007213 


22 






0.025031 


0.027608 


0.029788 


23 






-0.005219 


-0.006398 


-0.007213 


24 






0.013840 


0.012139 


0.012059 


25 






0.001828 


0.002663 


0.002865 


26 






0.025031 


0.027608 


0.029788 


27 






0.004266 


0.006421 


0.007184 


28 






0.000679 


0.004521 


0.006052 


29 






-0.004266 


-0.006421 


-0.007184 


30 






0.000679 


0.004521 


0.006052 
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